The Sensitivity of Multidimensional Nova Calculations to the Outer Boundary Conditions

Multidimensional reactive flow models of accreted hydrogen rich envelopes on top of degenerate cold white dwarfs are very effective tools for the study of critical, non spherically symmetric, behaviors during the early stages of nova outbursts. Such models can shed light both on the mechanism responsible for the heavy element enrichment observed to characterize nova envelope matter and on the role of perturbations during the early stages of ignition of the runaway. The complexity of convective reactive flow in multi-dimensions makes the computational model itself complex and sensitive to the details of the numerics. In this study, we demonstrate that the imposed outer boundary condition can have a dramatic effect on the solution. Several commonly used choices for the outer boundary conditions are examined. It is shown that the solutions obtained from Lagrangian simulations, where the envelope is allowed to expand and mass is being conserved, are consistent with spherically symmetric solutions. In Eulerian schemes which utilize an outer boundary condition of free outflow, the outburst can be artificially quenched.Comment: 12 Pages 3 figures; Accepted for publication in the Astrophysical Journa

Nucleosynthesis in novae

Peer ReviewedPostprint (published version

Modeling Core-Collapse Supernovae in 3-Dimensions

We present the first complete 3-dimensional simulations of the core-collapse of a massive star from the onset of collapse to the resultant supernova explosion. We compare the structure of the convective instabilities that occur in 3-dimensional models with those of past 2-dimensional simulations. Although the convective instabilities are clearly 3-dimensional in nature, we find that both the size-scale of the flows and the net enhancement to neutrino heating does not differ greatly between 2- and 3-dimensional models. The explosion energy, explosion timescale, and remnant mass does not differ by more than 10% between 2- and 3-dimensional simulations.Comment: 5 pages text, 3 separate figures (see http://qso.lanl.gov/~clf for more info), accepted by Ap

Gravity Waves in the Sun

We present numerical simulations of penetrative convection and gravity wave excitation in the Sun. Gravity waves are self-consistently generated by a convective zone overlying a radiative interior. We produce power spectra for gravity waves in the radiative region as well as estimates for the energy flux of gravity waves below the convection zone. We calculate a peak energy flux in waves below the convection zone to be three orders of magnitude smaller than previous estimates for m=1. The simulations show that the linear dispersion relation is a good approximation only deep below the convective-radiative boundary. Both low frequency propagating gravity waves as well as higher frequency standing modes are generated; although we find that convection does not continually drive the standing g-mode frequencies.Comment: 22 pages, 14 figures, submitted to MNRA

Calibrating the Mixing Length Parameter for a Red Giant Envelope

Two-dimensional hydrodynamical simulations were made to calibrate the mixing length parameter for modeling red giant's convective envelope. As was briefly reported in Asida & Tuchman (97), a comparison of simulations starting with models integrated with different values of the mixing length parameter, has been made. In this paper more results are presented, including tests of the spatial resolution and Large Eddy Simulation terms used by the numerical code. The consistent value of the mixing length parameter was found to be 1.4, for a red giant of mass 1.2 solar-mass, core mass of 0.96 solar-mass, luminosity of 200 solar-luminosities, and metallicity Z=0.001.Comment: 18 pages, 1 table, 13 figures. Accepted for publication in Ap.

Presolar Grains from Novae

We report the discovery of five SiC grains and one graphite grain isolated from the Murchison carbonaceous meteorite whose major-element isotopic compositions indicate an origin in nova explosions. The grains are characterized by low 12C/13C (4-9) and 14N/15N (5-20) ratios, large excesses in 30Si (30Si/28Si ratios range to 2.1 times solar) and high 26Al/27Al ratios. These isotopic signatures are theoretically predicted for the ejecta from ONe novae and cannot be matched by any other stellar sources. Previous studies of presolar grains from primitive meteorites have shown that the vast majority formed in red giant outflows and supernova ejecta. Although a classical nova origin was suggested for a few presolar graphite grains on the basis of 22Ne enrichments, this identification is somewhat ambiguous since it is based only on one trace element. Our present study presents the first evidence for nova grains on the basis of major element isotopic compositions of single grains. We also present the results of nucleosynthetic calculations of classical nova models and compare the predicted isotopic ratios with those of the grains. The comparison points toward massive ONe novae if the ejecta are mixed with material of close-to-solar composition.Comment: 20 pages, 5 figures, 1 table. ApJ, in pres

On Validating an Astrophysical Simulation Code

We present a case study of validating an astrophysical simulation code. Our study focuses on validating FLASH, a parallel, adaptive-mesh hydrodynamics code for studying the compressible, reactive flows found in many astrophysical environments. We describe the astrophysics problems of interest and the challenges associated with simulating these problems. We describe methodology and discuss solutions to difficulties encountered in verification and validation. We describe verification tests regularly administered to the code, present the results of new verification tests, and outline a method for testing general equations of state. We present the results of two validation tests in which we compared simulations to experimental data. The first is of a laser-driven shock propagating through a multi-layer target, a configuration subject to both Rayleigh-Taylor and Richtmyer-Meshkov instabilities. The second test is a classic Rayleigh-Taylor instability, where a heavy fluid is supported against the force of gravity by a light fluid. Our simulations of the multi-layer target experiments showed good agreement with the experimental results, but our simulations of the Rayleigh-Taylor instability did not agree well with the experimental results. We discuss our findings and present results of additional simulations undertaken to further investigate the Rayleigh-Taylor instability.Comment: 76 pages, 26 figures (3 color), Accepted for publication in the ApJ

A New 17F(p,gamma)18Ne Reaction Rate and Its Implications for Nova Nucleosynthesis

Proton capture by 17F plays an important role in the synthesis of nuclei in nova explosions. A revised rate for this reaction, based on a measurement of the 1H(17F,p)17F excitation function using a radioactive 17F beam at ORNL's Holifield Radioactive Ion Beam Facility, is used to calculate the nucleosynthesis in nova outbursts on the surfaces of 1.25 and 1.35 solar mass ONeMg white dwarfs and a 1.00 solar mass CO white dwarf. We find that the new 17F(p,gamma)18Ne reaction rate changes the abundances of some nuclides (e.g., 17O) synthesized in the hottest zones of an explosion on a 1.35 solar mass white dwarf by more than a factor of 10,000 compared to calculations using some previous estimates for this reaction rate, and by more than a factor of 3 when the entire exploding envelope is considered. In a 1.25 solar mass white dwarf nova explosion, this new rate changes the abundances of some nuclides synthesized in the hottest zones by more than a factor of 600, and by more than a factor of 2 when the entire exploding envelope is considered. Calculations for the 1.00 solar mass white dwarf nova show that this new rate changes the abundance of 18Ne by 21%, but has negligible effect on all other nuclides. Comparison of model predictions with observations is also discussed.Comment: 20 pages, 6 figures, accepted for publication in Ap

Nucleosynthesis in ONeMg Novae: Models versus Observations to Constrain the Masses of ONeMg White Dwarfs and Their Envelopes

Nucleosynthesis in ONeMg novae has been investigated with the wide ranges of three parameters, i.e., the white dwarf mass, the envelope mass at ignition, and the initial composition. A quasi-analytic one-zone approach is used with an up-to-date nuclear reaction network. The nucleosynthesis results show correlation with the peak temperatures or the cooling timescales during outbursts. Among the combinations of white dwarf and envelope masses which give the same peak temperature, the explosion is more violent for a lower white dwarf mass owing to its smaller gravitational potential. Comparison of the nucleosynthesis results with observations implies that at least two-third of the white dwarf masses for the observed ONeMg novae are $\simeq 1.1 M_\odot$, which are significantly lower than estimated by previous hydrodynamic studies but consistent with the observations of V1974 Cyg. Moreover, the envelope masses derived from the comparison are $\gtrsim 10^{-4} M_\odot$, which are in good agreement with the ejecta masses estimated from observations but significantly higher than in previous hydrodynamic studies. With such a low mass white dwarf and a high mass envelope, the nova can produce interesting amounts of $\gamma$-ray emitters $^7$Be, $^{22}$Na, and $^{26}$Al. We suggest that V1974 Cyg has produced $^{22}$Na as high as the upper limit derived from the COMPTEL survey. In addition, a non-negligible part of the Galactic $^{26}$Al may originate from ONeMg novae, if not the major contributors. Both the future INTEGRAL survey for these $\gamma$-ray emitters and abundance estimates derived from ultraviolet, optical, and near infrared spectroscopies will impose a severe constraint on the current nova models.Comment: 21 pages, 23 figures, to appear in the Astrophysical Journal, Vol. 523, No.1, September 20, 1999; preprint with embedded images can be obtained from http://th.nao.ac.jp/~wanajo/journal/onenova.p